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Abstract 

The field-theoretical renormalization group (RG) approach in three dimen- 
sions is used to estimate the universal critical values of renormalized coupling 
constants g§ and g& for the 0(n)-symmetric model. The RG series for g$ and 
gs are calculated in the four-loop and three-loop approximations, respectively, 
and then resummed by means of the Pade-Borel-Leroy technique. Under the 
optimal value of the shift parameter b providing the fastest convergence of the 
iteration procedure, numerical estimates for g§ are obtained with the accu- 
racy no worse than 0.3%. The RG expansion for g% demonstrates a stronger 
divergence, and results in considerably cruder numerical estimates. 
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I. INTRODUCTION 



The three-dimensional (3D) 0(n)-symmetric model plays a very important role in the 
theory of phase transitions. It describes critical phenomena in a variety of physical systems 
including Ising, XY-like and Heisenberg ferromagnets, simple fluids and binary mixtures, 
superconductors and Bose superffuids, etc. This model is also relevant to certain asymptotic 
regimes of the critical behavior of the quark-gluon plasma in quantum chromodynamics 
(n = 4) [0,^]. In the critical region, the n- vector model is known to be thermodynamically 
equivalent to the 3D Euclidean field theory of A<£> 4 type, and may be treated by the field- 
theoretical renormalization group (RG) technique which proved to be very efficient both for 
studying the qualitative features of phase transitions and calculating the critical exponents 
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On the other hand, for decades the influence of ordering fields upon the critical behavior 
of various systems attracted permanent attention, being of prime interest both for theorists 
and experimentalists. Recently, the free energy (effective action) and, in particular, higher- 
order renormalized coupling constants g<ik for the basic models of phase transitions became 
the target of intensive theoretical studies pH^. These constants are related to the non- 
linear susceptibilities X2k and enter the scaling equation of state, thus playing a key role at 
criticality. Along with critical exponents and critical amplitude ratios, they are universal, i.e. 
they possess, under T — > T c , numerical values that are not sensitive to the physical nature 
of the phase transition, depending only on the system dimensionality and the symmetry of 
the order parameter. 

Calculation of the universal critical values of g$, g$, etc. for the three-dimensional Ising 
model by a number of analytical and numerical methods showed that the field-theoretical 
RG approach in fixed dimensions yields the most accurate numerical estimates for these 
quantities. It is a consequence of a rapid convergence of the iteration schemes originating 
from renormalized perturbation theory. Indeed, the resummation of four- and five-loop RG 
expansions by means of the Borel-transformation-based procedures gave the values for 



which differ from each other by less than 0.5% [TgJT^], while the use of a resummed three 



loop RG expansion enabled one to achieve an apparent accuracy no worse than 1.6% [|7|,|T7 



Moreover, the field-theoretical RG approach turns out to be powerful enough even in two 
dimensions: properly resummed four-loop RG expansions lead to fair numerical estimates 
for the critical exponents [0] and the renormalized coupling constant gl P^J of a 2D Ising 



model, and give reasonable results for its random counterpart [25,26]. It is natural, therefore, 



to use the field theory for a calculation of renormalized higher-order coupling constants for 
the 3D n-vector model. In this paper, the 3D RG expansion for the renormalized coupling 
constants g§ and g$ will be calculated, and the numerical estimates for their universal critical 
values will be obtained. 



II. RG EXPANSIONS FOR THE SEXTIC AND OCTIC COUPLING CONSTANTS 

Within the field-theoretical language, the 3D 0(n)-symmetric model in the critical region 
is described by Euclidean scalar field theory with the Hamiltonian 



H 



d 6 x 



2 \2 



(1) 



where a bare mass squared is proportional to T — Tjp\ Tj°) being the phase transition 
temperature in the absence of the order parameter fluctuations. Taking fluctuations into 
account results in renormalizations of the mass m — > m, the field tp — > ip^, and the coupling 
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constant A — > mg^. Moreover, thermal fluctuations give rise to many-point correlations 
< ty(x\)ip{%2)...ip{x2k) > and, correspondingly, to higher-order terms in the expansion of the 
free energy in powers of the magnetization M: 



F(M, m) = F(0, m) + J2 r 2fc^/ 

fc=i 



2k 



(2) 



In the critical region, the coefficients T 2 k, being one-particle irreducible 2fc-point vertices 
taken at zero external momenta, demonstrate well-known scaling behavior: 



2k — 92k,™ , 



(3) 



where n is a Fisher exponent, and g 2 h are some constants. Let us set as usually g 2 = 
1/2. Then g A , g & , g%,... will acquire universal values. The asymptotic critical values of g±, 
gl(n), determine the critical exponents and other universal quantities, thus playing a very 
important role in the theory. The numbers gl(n) have been found by resummation of the 
six-loop expansion for the RG /3-function E ,f||§[7| , from strong-coupling series [^7]] , by lattice 
calculations [2TJ, from the e-expansion |28], and are known today with an accuracy which 
may be considered rather high. 

The universal values of higher coupling constants g$, g£, etc. determine the structure of 
the free energy F(M, m) under strong critical fluctuations. In fact, the Taylor expansion of 
the scaling function contains the ratios g 2k / / (gl) k ~ 1 , which may be easily shown by replace- 
ment of the magnetization M in Eq.(|2]) by the dimensionless variable z = M^Jg^/m 1 ^: 



ev \ i7 /n \ m 3 I z 2 4 g 6 6 g 8 g \ 
F(z, m) - F(0, m) = — — + z + -zz° + -^z° + ... . 

9i V 2 gi g% J 



(4) 



Moreover, via g 2 k, the non-linear susceptibilities \2k can be expressed. For \i an d Xe 
corresponding formulas are as follows: 



d 3 M 



Xe 



H=0 



d 5 M 
dH 5 



720 X lm- 6 (8gl-g 6 ). (5) 



H=0 



Their inversion gives the relations 



9i 



m z XA 



96 



m 6 (10x4 



X6X2) 



720x1 



(6) 



which are widely used for extraction of numerical values of renormalized coupling constants 
from the results of lattice calculations [14,16,21,22,29,30]. 

The method of calculating the RG series for the g§ and g$ we use here is straightfor- 
ward. Since in three dimensions higher-order bare couplings are irrelevant in RG sense, the 
renormalized perturbative series to be found can be obtained from conventional Feynman 
graph expansions for the six-point and eight-point vertices in terms of the only bare coupling 
constant - A. In the course of calculations the tensor structure of these vertices, 



= y^(<W<$ 7 <5(W + 14 transpositions^ 6 . 



(7) 



105 



(SapdjgSnvSpp + 104 transpositions)! g 



(8) 
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should be taken into account. In its turn, A may be expressed perturbatively as a function 
of the renormalized coupling constant g^. Substituting corresponding power series for A into 
original expansions, we can obtain the RG series for g 6 and g$. The one-, two-, three- and 
four- loop contributions to g§ are formed by one, three, 16, and 94 one-particle irreducible 
Feynman graphs, respectively. Their calculation gives: 



9/AZ 2 



3 r 



96 = - 



tt \ m 



n + 26 9 n 2 + 340 n + 2324 f XZ 2 



27 



1627T 



\ m 



+ (0.0056289546468 n 3 + 0.28932672886 n 2 + 4.0404241235 n + 16.204286853) 



(XZ 2 \ 2 
\ m J 



-(0.001493126 n 4 + 0.09961447 n 3 + 2.152320 n 2 + 18.330704 n + 52.830284) 



'xz 



2\ 3- 



m 



(9) 



The perturbative expansion for A emerges directly from the normalizing condition A = 
mZ^Z~ 2 g^ and the known series for Z 4 ||: 

1 n + 8 3 n 2 + 38 n + 148 2 

Z4 = 1 + ^T^ 4 + 12^ 9 * 

+ (0.0040314418 n 3 + 0.0679416657 n 2 + 0.466356233 n + 1.240338484)# 3 . (10) 
Combining these expressions, we obtain 
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ge = -gl 

71 



n + 26 17 n + 226 

+ (0.000999164 n 2 + 0.14768927 n + 1.24127452)^ 



27 



817T 



- (-0.00000949 n 3 + 0.00783129 n 2 + 0.34565683 n + 2.14825455)^' 



(11) 



In the case of g$, the one-, two-, and three-loop contributions are given by one, five, 
and 36 Feynman graphs, respectively. Corresponding "bare" and renormalized perturbative 
expansions are found to be: 



9s 



81 (\Z 
2tt \ m 



2\ 4 



n + 80 405 n 2 + 35626 n + 342320 XZ 2 



13122vr 



ni 



+ (0.0046907955 n 3 + 0.463650683 n 2 + 8.86811653 n + 45.4769028) 



'XZ 2 



(12) 



81 4 



n + 80 81 n 2 + 7114 n + 134960 



81 



13122tt 



-94 



+ (0.00943497 n 2 + 0.60941312 n + 7. 15615323) g\ 



(13) 



In Sec. Ill, the series Eqs.flTTD, (|13|) will be used for estimation of the universal numbers 
9t and 9l 
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III. RESUMMATION AND NUMERICAL ESTIMATES 



Being a field-theoretical perturbative expansions the series of equations (ITlf) , (|13|) have 
factorially growing coefficients, i. e., they are divergent (asymptotic). Hence, direct sub- 
stitution of the fixed point value g\ into them would not lead to satisfactory results. To 
get reasonable numerical estimates for g$ and g$, some procedure making these expansions 
convergent should be applied. As is well known, the Borel-Leroy transformation 



OO ^ CO 

f(x) = £ ca l = / *e*F(xt)dt, F(y) = £ j-^rj, (14) 

diminishing the coefficients by the factor (i + b)\, can play a role of such a procedure. Since 
the RG series considered turns out to be alternating the analytical continuation of the Borel 
transform may be then performed by using Pade approximants. 

Let us discuss first the estimation of the sextic coupling constant g§. With the four-loop 
expansion ([TT|) in hand, we can construct, in principle, three different Pade approximants: 
[2/1], [1/2], and [0/3]. To obtain proper approximation schemes, however, only diagonal 
[L/L] and near-diagonal Pade approximants should be employed |31]| . That is why, further, 
when estimating we limit ourselves with approximants [2/1] and [1/2]. Moreover, the 
diagonal Pade approximant [1/1] will be also dealt with, although this corresponds, in fact, 
to the usage of the lower-order, three-loop RG approximation. 

The algorithm of estimating g§ we use here is as follows. Since the Taylor expansion for 
the free energy contains as coefficients the ratios i?2fc = fi^fc/fi^ -1 rather than the renormalized 
coupling constants themselves, we work with the RG series for R§. It is resummed in 
three different ways based on the Borel-Leroy transformation and the Pade approximants 
just mentioned. The Borel-Leroy integral is evaluated as a function of the parameter h 
under g^ = g\. For the fixed point coordinate gl, the values given by the resummed six- 
loop RG expansion for the /3-function are adopted |||7|, which are believed to be the most 
accurate estimates available today. The optimal value of b providing the fastest convergence 
of the iteration scheme is then determined. It is deduced from the condition that the Pade 
approximants employed should give, for b = b opt , the values of i?g which are as close as 
possible to each other. Finally, the average over three estimates for i?g is found and claimed 
to be a numerical value of this universal ratio. 

To obtain an idea about how such a procedure works, let us use Table I, where the 
results of corresponding calculations for n = 1, 3, and 10 are presented. It is seen that 
for n = 1 and 3, b opt , providing a coincidence of the estimates given by all three working 
Pade approximants, is equal to 1.24. For n = 10, b opt , fixed by the approximants [1/1] 
and [2/1], is equal to 1, whereas the third approximant ([1/2]) at b = 1 is spoiled by a 
positive axis pole. Nevertheless, the numerical estimate given by this approximant under 
the nearest "safe" (integer) value of b {b = 2) turns out to be very close to that predicted 
by the pole free approximants for b opt . Moreover, as is seen from Table I, with increasing 
n numerical estimates for g$ become less dependent on b, i. e., their sensitivity to the type 



of resummation decreases. This is not surprising. The point is that the RG expansion (|TT|) 
becomes less divergent when n grows up. To make this property obvious, let us replace g^ 
in Eq.([y]) by the effective coupling constant 

9 = ~^-94, (15) 

that is known to be only weakly dependent on n: it varies from 1.415 to 1 when n goes from 
1 to infinity Then we obtain 
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_ 8vr 2 (n + 26) 3 
96 ~ 3(n + 8) 3 9 



2(17n + 226) 
3(n + 8)(n + 26) 
-0.0638n 3 



1.065025™ 2 + 157.42454™ + 1323.09596 



(n + 8) 2 (n + 26) 
52.4510n 2 + 2314.9897n + 14387.6460 



n + 



in 



26) 



g" 



(16) 



One can see now that all the terms in the RG expansion for g 6 (in square brackets), apart 
from the first one, decrease monotonically when n — > oo. This implies that the larger n is 
the smaller the contribution of the higher-order terms and, correspondingly, the better the 
approximating properties of this series. 

This conclusion is definitely confirmed by Table II. It contains numerical estimates for 
#g resulting from the four-loop RG expansion resummed by the Pade-Borel-Leroy technique 
described above (column 3) and their analogs given by the Pade-Borel resummed three-loop 
RG series M (column 4). As is seen, with increasing n the difference between the four- loop 
and three-loop estimates rapidly diminishes. Being small (0.9 %) even for n — 1, it becomes 
negligible at n — 10 and practically disappears for n > 14. 

How close to the exact values of may the numbers in column 3 be? To clear up this 
point, let us compare our four-loop estimate for R§ at n = 1 with those obtained recently by 
an analysis of the five-loop scaling equation of state for the 3D Ising model ]T9],|52]]. Guida 



and Zinn- Justin obtained i?g = 1.644 and, taking into account some additional information, 
Rq = 1.643, while our estimate is R$ = 1.648. Keeping in mind that the exact value of R§ 
should lie between the four-loop and five- loop estimates (the RG series is alternating), our 
estimate can differ from the exact number by no more than 0.3 %. Since for n > 1 the RG 
expansion (|ll|) was argued to provide the better numerical estimates than in the Ising case, 
this value (0.3 %) may be referred to as an upper bound for the deviation of the numbers 
in column 3 of Table II from their exact counterparts. 

It is interesting to compare our estimates for with those obtained by other methods. 
Since 1994, the universal values of the sextic coupling constant for the 3D 0(n)-symmetric 
model were estimated by solving the exact RG equations [11], by lattice calculations [14 . 
and by a constrained analysis of the e-expansion [p3|1 ; corresponding results are collected in 
columns 5, 6, and 7 of Table II, respectively. As is seen, they are, in general, in accord with 
ours. 

A less optimistic situation takes place in the case of the octic coupling constant g$. The 
RG expansion Eq . (|13|) is shorter than Eq . ([TT|) , and stronger diverges. Moreover, the second 
term in this series, along with the first one, remains finite under n — > oo. It becomes obvious 
if one replaces g^ by g: 



3vr 3 (n + 80) 4 

- g 



in + 



:| 81 n 2 + 7114 n + 134960 

81(n + 80)(n + 8) 9 



30.1707 n 2 + 1948.7519 n + 22883.6021 2 
(n + 80)(n + 8) 2 9 



(17) 



In addition, the RG series for g 8 has an unusual feature: when n — > oo, the first and 
second terms tend to compensate each other, making their mutual contribution small and 
increasing the role of the higher-order terms. That is why numerical estimates resulting 
from the expansion ( |TTD are expected to be substantially cruder than those given by the 
series ( |i~6|) both for small and large values of n. 

In order to estimate gl{n), we resum the RG expansion for g 8 by the Pade-Borel-Leroy 
technique using the diagonal Pade approximant [1/1]. Other Pade approximants, [0/2] and 
[0/1], are ignored, since they turn out to lead to quite unsatisfactory numerical results. 
Dealing with a single Pade approximant, in some condition we need to fix the optimal value 
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of the shift parameter b. For the three-dimensional Ising model the estimate g$ = 0.825 was 
recently found |19[]. This number has been extracted from the five-loop RG expansion, and 
may be considered the most accurate known up to the present. It is natural therefore to 
tune, by proper choice of b, a numerical value of g|(l) gi ven by the resummed three-loop 
RG series with the best estimate available. Such a procedure leads to b opt = 40, and this 
number is adopted as optimal in the course of evaluation of for arbitrary n. 

The results of our calculations are collected in Table III, where the estimates for gg{n), 
obtained by a constrained analysis of the e-expansion | 2~3fl , by approximate solution of the 
exact RG equations |TT|, and given by the 1/n-expansion technique, are also presented for 



comparison. As seen, for n > 8 the numbers originating from two field-theoretical approaches 
- g-expansion in three dimensions and e-expansion - agree quite well. However, for smaller 
n, especially for n = 2, differences between them turn out to be rather large. This is 
not surprising since overly short perturbative expansions for g 8 are available both in 3 and 
4 — e dimensions and they demonstrate a strong divergence preventing accurate numerical 
estimates from being obtained. At the same time, our three-loop RG estimates are believed 
to be closer to the true critical values of g$ than those given by the e-expansion, because in 
three dimensions we have longer perturbative series. A fair agreement between our results 
and the numbers emerging from the exact RG equations (see Table III) may be considered 
as an argument in favor of this belief. 



IV. CONCLUSION 



To summarize, we have calculated the RG expansions for renormalized coupling con- 
stants #6 and g$ of the 3D n-vector model in four-loop and three-loop orders, respectively. 
Resummation of the RG series by the Pade-Borel-Leroy method has enabled us to obtain nu- 
merical estimates for the universal critical values of these quantities for arbitrary n. Having 
analyzed the sensitivity of the RG estimates for g§ to the type of resummation procedure 
and a character of their dependence on the order of the RG approximation, an apparent 
accuracy of these numbers has been argued to be no worse than 0.3%. Numerical estimates 
for gl turned out to be much less accurate because of a smaller length and stronger diver- 
gence of the RG expansion obtained. They were found to be consistent, in general, with the 
values of g§ deduced from the exact RG equations and, for n > 8, with those given by a 
constrained analysis of corresponding e-expansion. 
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TABLES 



TABLE I. The values of #g for n = 1, 3, and 10 obtained by means of the Pade-Borel-Leroy 
technique for various b within three- loop (approximant [1/1]) and four-loop (approximants [1/2] 
and [2/1]) RG approximations. The estimates for several values of b in the middle lines are absent 
because corresponding Pade approximant turns out to be spoiled by a positive axis pole. 



b 





1 


1.24 


2 


3 


4 


5 


7 


n = 1 


[1/1] 


1.576 


1.604 


1.6089 


1.621 


1.633 


1.641 


1.648 


1.656 


[1/2] 






1.6084 


1.600 


1.595 


1.592 


1.590 


1.587 


[2/1] 


1.639 


1.613 


1.6084 


1.596 


1.583 


1.573 


1.566 


1.555 


n = 3 


[1/1] 


0.937 


0.949 


0.95133 


0.957 


0.962 


0.966 


0.969 


0.973 


[1/2] 






0.95133 


0.948 


0.946 


0.944 


0.944 


0.942 


[2/1] 


0.964 


0.953 


0.95133 


0.946 


0.941 


0.937 


0.934 


0.930 


n = 10 


[1/1] 


0.2338 


0.23515 




0.2360 


0.2366 


0.2370 


0.2373 


0.2377 


[1/2] 








0.2348 


0.2346 


0.2345 


0.2344 


0.2342 


[2/1] 


0.2359 


0.23515 




0.2346 


0.2342 


0.2339 


0.2337 


0.2334 
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TABLE II. Our estimates of universal critical values of the renormalized sextic coupling con- 
stant for the 3D n- vector model (column 3). The fixed point coordinates g* are taken from Ref.[3] 
(1 < n < 3) and Ref.[7] (4 < n < 40). The g§ estimates extracted earlier from Pade-Borel 
resummed three- loop RG expansion (column 4), from the exact RG equations (column 5), ob- 
tained by the lattice calculations (column 6), and resulting from a constrained analysis of the 
e-expansions (column 7) are presented for comparison. Column 8 contains the values of g§ given 
by the 1/n-expansion technique. 



n 


n* 

9 


n* 
96 


n* 71 

9e lill 


n* [1 1 1 
96 lliill 


1141 
96 v3 


96 l|£J 


96 \ l / n ) 




9 

z 


O 


/I 


K 



D 


7 


Q 
O 


i 
i 


1 41 ^ 


i fins 

l.UUO 


1 R99 


1 ^9 
i.oz 


i.yzi z*± i 


i.uuy^y j 




9 


1 4f")fi 


1 998 


l.ZOO 


114 


1 97f9 c i N l 
l.Z / ^ZO J 


1 91 ^7^ 




Q 
O 


1 *}Q9 

i.oyz 


u.yoi 


u.yoo 


U.oo 


u.yo^zu j 


u.yoi ^4to j 




A 

-± 


1.0 / 4LO 


n 747 


n 7M 

U. / 01 


U.uo 


n Ron K\ 

U.OZ^IO J 


n 79^90"! 

u. / zo^zy j 


1 fi44Q 


K 
o 


1 .OJUO 


U.OyU 


u.oyy 








1 0^98 

1 .UOZO 


u 


l.OOOU 




n 48^ 








n 7^1 1 

u. / Oil 


7 


l.OZl 


u.oyo 


u.oyo 








n ^71 
u.oo * i 


o 
o 


1.0U4LO 


n ^9Q 


n ^1 

U.OOl 






n ^1 0^4^ 
u.oiy ^ j 


n 41 1 9 

U.4:l 1Z 


9 


1.289 


0.277 


0.278 








0.3249 


10 


1.2745 


0.235 


0.236 








0.2632 


12 


1.2487 


0.174 


0.175 








0.1828 


14 


1.2266 


0.134 


0.134 








0.1343 


16 


1.2077 


0.105 


0.105 






0.1032(4) 


0.1028 


18 


1.1914 


0.0845 


0.0847 








0.0812 


20 


1.1773 


0.0693 


0.0694 








0.0658 


24 


1.1542 


0.0487 


0.0488 








0.0457 


28 


1.1361 


0.0360 


0.0361 








0.0336 


32 


1.1218 


0.0276 


0.0276 






0.0275(1) 


0.0257 


36 


1.1099 


0.0218 


0.0218 








0.0203 


40 


1.1003 


0.0176 


0.0176 








0.0164 
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TABLE III. Three-loop RG estimates of universal critical values of the renormalized octic 
coupling constant g$ (column 2). The g$ estimates resulting from a constrained analysis of the 
e-expansion (column 3), from the exact RG equations (column 4) and given by the 1/n-expansion 
technique (column 5) are presented for comparison. 



n 


9s 


9s JfJ 


(ill 

9s \ u \ 


9s \ l / n ) 




z 


Q 
O 


4 


r 
O 


1 


U.ozo 


n on/n\ 

u.oz^y j 


n vol 
U. / zl 




Z 


0.000 


0.o3(31J 


0.343 




3 


0.168 


0.36(17) 


0.145 




4 


0.057 


0.15(13) 


0.042 


-2.151 


6 


-0.021 






-0.834 


8 


-0.034 


-0.03(2) 




-0.0388 


16 


-0.014 


-0.015(2) 




-0.0456 


32 


-0.0023 


-0.0023(1) 




-0.00395 


48 


-0.00062 


-0.00061(2) 




-0.00087 


64 


-0.00023 






-0.00029 


100 


-0.000046 


-0.000044(2) 


-0.000049 


-0.000052 
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